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ABSTRACT 

In classical invariant theory, the Grobner base of the ideal 
of syzygies and the normal forms of polynomials of invari- 
ants are two core contents. To improve the performance of 
invariant theory in symbolic computing of classical geome- 
try, advanced invariants are introduced via Clifford product 
[5] . This paper addresses and solves the two key problems in 
advanced invariant theory: the Grobner base of the ideal of 
syzygies among advanced invariants, and the normal forms 
of polynomials of advanced invariants. These results beauti- 
fully extend the straightening of Young tableaux to advanced 
invariants. 
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1. INTRODUCTION 

In traditional analytic approach to classical geometry, coor- 
dinates are introduced to represent points in the geometric 
space, and equations of the coordinates are used to define 
constraints among the points, forming a representation of 
higher dimensional objects such as curves, surfaces, etc. Ba- 
sic manipulations of coordinates include addition and multi- 
plication, resulting in polynomials in the coordinates. Since 
the coordinates of generic points are independent, and the 
multiplication of coordinate variables are commutative, nor- 
malization of the polynomials in the coordinates is very easy. 
The normal forms of the polynomials are required in many 
manipulations, e.g., division among polynomials. 

Another analytic approach to classical geometry, dating back 
to Euclid, is to use geometric invariants such as lengths, an- 
gles, areas, etc. A typically algebraic system of geometric 
invariants is a polynomial ring generated by basic invariants. 
In such a system, a vector variable in a linear space is used 
to represent a point or direction in classical geometry, the 
inner product of a vector with itself represents the squared 



length of the vector, the inner product of two unit vectors 
represents the cosine of the angle between them, etc. Such 
operators among vectors generate a set of basic invariants, 
and the polynomials in these basic invariants are advanced 
invariants. 

Although the multiplication of invariants are commutative, 
the basic invariants generated by generic vector variables of 
the linear space are not independent, and there are poly- 
nomial relations among them, called syzygy relations. The 
dependency is largely caused by the dimension constraint 
of the linear space upon vectors. While the dimension con- 
straint can be easily reflected by the number of coordinates 
introduced to represent a point and the independency among 
the coordinates, for basic invariants generated by the points, 
fully representing the dimension constraint is by no means 
trivial. Classical invariant theory studies the generators of 
invariants, the syzygy relations among the basic invariants, 
and the normal forms of advanced invariants as polynomials 
in the basic ones [11] . [12] , 

In symbolic geometric computing, both the coordinate ap- 
proach and the basic invariant approach encounter the dif- 
ficulty of very big polynomial size, in particular in the mid- 
dle of symbolic manipulations. In ,.5j, a recipe to allevi- 
ate the difficulty is proposed, called long geometric product, 
BREEFS, and Clifford factorization, among which the long 
geometric product (or Clifford product) is the foundation. 
The idea is to convert polynomials of basic invariants into 
advanced invariants, converse to the approach of classical 
invariant theory, by means of an associative and multilin- 
ear product among the vector variables representing points. 
The associativity of the product and the symmetries within 
a long bracket provide powerful manipulations that cannot 
be done with basic invariants, nor with coordinates. This is 
a top-down approach to advanced invariants [6], while the 
classical invariant theory is a bottom-up approach. 

Dealing with the syzygy relations among advanced invari- 
ants and finding the normal forms of polynomials in ad- 
vanced invariants are two fundamental tasks in such "ad- 
vanced invariant theory". The 2D case is easy, while higher 
dimensional cases are difficult. Little advance has been 
achieved in six years since the publication of [5] in 1997. 

In this paper, the two fundamental problems are solved for 
the advanced invariant theory of 3D orthogonal geometry: 
the Grobner base of the syzygy ideal of "long brackets", and 



the normal forms of Clifford bracket polynomials, ft turns 
out that the normal forms of such bracket polynomials are 
surprisingly "beautiful". The description is the following. 

In classical invariant theory for (n— 1)D projective geometry, 
the basic invariants are brackets of length n, or in coordinate 
form, the n x n determinants formed by the homogeneous 
coordinates of n vector variables. A bracket polynomial is in 
normal form if when each term is up to coefficient written in 
Young tableau form, the entries in each row are increasing, 
while the entries in each column are non-decreasing [14] . For 
example for vector variables vi -< V2 -< . . - -< v m , a bracket 
monomial [v ill v H2 • ■ • v ilc ][vi 21 Vi 22 ■ • ■ v; 2c ] ■ • ■ [w irl w ir2 ■ ■ ■ 
v, rc ], where the Vj. fe are repetitive selections of the m vector 
variables, is normal if and only if in 



Vin v ila 

V »21 V »22 



(1.1) 



v, rl v ir2 ■ • ■ v irc 

Viji X v i 3 2 X < v bc while Vi u X v i2k X • . . X v irk . 

In the advanced invariant theory for 3D orthogonal geom- 
etry, each "elementary" advanced invariant is a bracket of 
length > 1, whose entries are vector variables representing 
3D points. In a bracket monomial, different brackets many 
have different lengths, and a bracket monomial is in normal 
form if and only if not only the entries in each row are in- 
creasing, the entries in each column are non-decreasing, but 
all the entries in the tableau after removing the first column, 
are non-decreasing. For example if (|1.1[) is normal in this set- 
ting, then the sequence v in v il2 • • • v ilc Vi 21 v i22 • • • Vi 2e • • • 
Vi rl Vi r2 ■ ■ • Vi rc is non-decreasing, where Vi fc denotes that 
Vi k does not occur in the sequence. 



This paper is organized as follows. Section 2 introduces 
orthogonal geometric invariants by quaternions, Clifford al- 
gebra and bracket algebra. Section 3 introduces the main 
results in 7: on vector-variable polynomials and basics of 
"advanced bracket algebra". Section 4 provides the Grobner 
base and normal forms of long brackets in the multilinear 
case. Section 6 extends the results to general case, by means 
of the square-free vector-variable polynomials introduced in 
Section 5. Section 7 proposes a normalization algorithm for 
bracket polynomials. 

2. QUATERNIONS, CLIFFORD ALGEBRA, 
AND BRACKET ALGEBRA 

In the vector algebra over R 3 , there are three multilinear 
products among vectors: (i) the inner product of two vec- 
tors, (ii) the cross product of two vectors, (iii) the hybrid 
product of three vectors. None of them can be extended to 
include more vectors while preserving the associativity. 

The quaternionic product, on the other hand, is associative 
while still multilinear. Let q represent the quaternionic con- 
jugate of quaternion q. Among quaternions, a vector v refers 
to a pure imaginary quaternion, i.e., v = — v, and a scalar 
s refers to a real quaternion, i.e., s — s. All vectors span a 
3D real inner-product space with metric diag(— 1, — 1, — 1), 
denoted by R~ 3 . 



We always use juxtaposition of elements to denote their 
quaternionic product. The inner product of two vectors 
v« , Vj is defined by 

[viVi] := (Vivj + VjV()/2. (2.1) 

The cross product of two vectors Vj, Vj is defined by 

v 4 x Vj := (v*Vj - VjVj)/2. (2.2) 

The result is a vector, so its inner product with a third 
vector Vfc is a scalar. Define the hybrid product of three 
vectors v< , v,- , Vfe by 

[v»VjVfc] := (viVjVfc - VfeVjVi)/2. (2.3) 

Then [v 4 VjV fc ] = [(vj x v J )v fc ]. 

The vector algebra over R -3 is included in the quaternions. 
The latter is equipped with a powerful product, the quater- 
nionic product, making it possible to use quaternions to rep- 
resent 3D orthogonal transformations pQ. 

The magnitude of a quaternion q is y'qq. A quaternion q 
is said to be unit if qq = 1. Let q be a unit quaternion, 
and v be a vector. The conjugate adjoint action of q on v 
is defined by 



Ad q (v) := qvq. 



(2.4) 



Since L4dq(vi)Ad q (v2)] = [V1V2] for any two vectors vi, V2, 
Adq realizes an orthogonal transformation in R~ 3 . A classi- 
cal result states that in fact all orthogonal transformations 
in R -3 are realized in this way, and two different unit quater- 
nions realize the same orthogonal transformation if and only 
if they differ by sign. 

For a quaternion Q, the bracket [Q] is its scalar part: 

[Q]:=(Q + Q)/2. (2.5) 

The axis of Q is the vector part of Q: 

A(Q) := (Q - Q)/2. (2.6) 

In particular, A(viV2) = vi x V2. We interpret them in 
geometrical terms below. 

For a unit vector vi, Ad vi realizes the reflection with re- 
spect to the plane normal to vi . In general, for unit vectors 
vi, V2, . . . , V2fc+i, Ad vlV2 ... V2k+1 realizes the reflection with 
respect to the plane normal to axis A(viV2 • ■ • V2fe+i), if the 
latter is nonzero. 

For two unit vectors vi,V2 that are linearly independent, 
Ad vlV2 realizes the rotation about the axis vi x V2: in the 
plane spanned by vi,V2, the rotation is from vi to the re- 
flection of vi with respect to V2, i.e., the angle of rotation is 
9 — 2Z(vi,V2). Furthermore, [V1V2] = cos(#/2). When we 
say "rotation V1V2", we mean the one induced by Ad v±V2 . 



In general, for unit vectors Vi, V2, . . . , V2fc, Ad vi ^ 



real- 



izes the rotation about the axis y4(viV2 • • • V2fc), if the latter 
is nonzero. The rotation is the composition of k rotations 
V1V2, V3V4, . . . , V2fc_iV2fc. If the angle of rotation is 9, then 



vi 



•v afc ] =cos(0/2), \A(vi---v 2k )\ = |sin(0/2)|. (2.7) 



Let A(viv 2 ■ • ■ v 2fc ) / 0. By [viV 2 ■ ■ ■ v 2fe+ i] = [A(viv 2 
■ • • v 2 fc)v2fc+i], we get 

[viv 2 • • ■ v 2fe+ i] = cosZ(A(viv 2 ---v 2fc ),v 2fe+ i) sin(0/2), 

(2.8) 

where 6 is the angle of rotation viv 2 • ■ ■ v 2 fe. In particular 
when k — 1, for linearly independent unit vectors vi,v 2 , 
sin(#/2) = sin Z(vi, v 2 ) equals the area of the parallelogram 
spanned by vi, v 2 , and cosZ(vi x v 2 ,V3) equals the height 
from the end of unit vector V3 to the plane spanned by 
vi,v 2 , so [V1V2V3] equals the volume of the parallelepiped 
spanned by vi,v 2 ,V3. 

In classical invariant theory, an algebraic invariant is a poly- 
nomial whose variables are basic invariants. In 3D orthogo- 
nal geometry, there are two kinds of basic invariants: [v<Vj] 
and [viVjVfc] for all vector variables Vj, v 3 *, Vfe. Given n vec- 
tor variables Vi, . . . , v n , the brackets Vj 2 ■ • ■ v Jm ] for ar- 
bitrary 1 < m < 00 and arbitrary repetitive selection of 
elements Vj 1 , v J2 , . . . , v Jm from the n variables, form an in- 
finite set of advanced algebraic invariants. That each of 
them is a polynomial of the [v;Vj] and [v;VjVfe] is guaran- 
teed by the following Caianiello expansion formulas [2], [6]: 
let V fe = v^Vjj ■ • ■ v Jfc , then 

[V 2 i] = E?i 2 (- !) 1 [vji Vj 4 ] [vj 2 v J3 • • • Vjv • • • v J2l } ■ 
A(V 2 i-i) = E( 2 ;- 2 ,i)hv 2 ,_J V 2i+i(i)] v 2;-i( 2 ); 

(2;-2,2)hV 2i 1 V 2!(1)]^4(V 2! (2)); 
[Val+l] = E(2;-2,3)hV 2i + 1 [ V 2* + l(l)][ V 2i+l(2)]> 

(2.9) 

where (i) (h, m— ft) h V m is a bipartition of the m elements 
in the sequence V m into two subsequences V m m and V m ( 2 ) 
of length /i and m — h respectively; (ii) in [V m m], the prod- 
uct of the h elements in the subsequence is denoted by the 
same symbol V m (i); (iii) the summation E(h m-fc)hv ' s 
over all such bipartitions of V m , and the sign of permu- 
tation of the new sequence V m m, V m (2) is assumed to be 
carried by the first factor [V m (i)] of the addend. 

While quaternions are sufficient for describing orthogonal 
transformations in 3D, they cannot be generalized to higher 
dimensions directly. In quaternions, the hybrid product 
[v;VjVfc] is a scalar. To make high-dimensional generaliza- 
tion this requirement must be removed, at the same time 
the property that this element be in the center of the alge- 
bra needs to be preserved. If we denote the quaternions by 
Q, then the above revision leads to a new algebra Q © lQ 
of dimension 8, where 1 := [V1V2V3] for three fixed vector 
variables that are linearly independent. This algebra is the 
Clifford algebra over R~ 3 . 

The formal definition of the Clifford algebra Cl(V n ) over an 
n-dimensional K-linear space V™ , where the characteristic of 
K is 7^ 2, is the quotient of the tensor algebra (^) V n over the 
ideal generated by elements of the form v ® v — Q(v) where 
Q is a K-quadratic form. The product induced from the 
tensor product is called the Clifford product, also denoted 
by juxtaposition of elements [3], [8]. 

When V n = R -3 , the quaternionic product of vectors is the 
image of their Clifford product under the homomorphism in- 
duced by mapping ttoa nonzero scalar. In Clifford algebra, 
1 is not a scalar, but called a pseudoscalar because it not 



only commutes with everything, but spans a ID real space 
containing all hybrid products. The concept quaternionic 
conjugate is replaced by the Clifford conjugate, which is the 
linear extension of the following operation: for any vectors 
Vi! , v j2 , . . . , v Jk , let V fe = v jl v j2 ■ ■ ■ v jk , then 

VT— (-l) fe Vt, where := v Ifc ■ ■ • v 12 v h , g 1Q s 

is the reversion of Vfe. 



With the Clifford conjugate, we can define the magnitude of 
Vfe, the conjugate adjoint action Adv k , the bracket [Vfe], the 
axis A(Vfc), together with the concepts not involving conju- 
gate: the inner product [vjVj], the cross product Vj x v, and 
the hybrid product [viVj-v/t], just the same as in the case of 
quaternions. The only difference is that since [viV^Vfe] is now 
a pseudoscalar, while j4(V2z-i) remains a vector, j4(V 2 ;) is 
not, but a pseudovector. Geometrically, when A(viv 2 ) 7^ 0, 
it represents the plane spanned by vectors vi, v 2 , or equiv- 
alently, the invariant plane of rotation Ad vlV2 supporting 

Vl, V2. 

We see that unlike the quaternions where there are only two 
kinds of objects of different dimensions: scalars which are 
usually called 0-D objects, and vectors which represent 1-D 
directions and so are usually called 1-D objects, in Clifford 
algebra C1(R -3 ) there are four kinds of objects of differ- 
ent dimensions. Besides scalars and vectors, there are pseu- 
dovectors which represent 2-D directions (planes) , and pseu- 
doscalars which represent 3-D orientations (spaces). This 
is the reason why C1(R~ 3 ) can be extended to higher di- 
mensions by being capable of discerning objects of different 
dimensions. 

To represent algebraic invariants in 3D orthogonal geome- 
try, using the quaternionic product or the Clifford product 
in the brackets does not make any difference. The geomet- 
ric interpretations (|2.7[l . (|2.8p and the Caianiello expansion 
formulas are identical for both products. This justifies the 
use of juxtaposition of elements to represent both products. 

The two kinds of basic invariants [viVj] and [viVjVfe] form a 
commutative ring, called inner-product bracket algebra. For- 
mally, given a set of n symbols M = {vi , . . . , v„}, two kinds 
of new symbols can be defined as following: (1) all 2-tuples 
of symbols selected repetitively from M, by requiring that 
each 2-tuple be symmetric with respect to its two entries; 
such a 2-tuple is denoted by [v;Vj]. (2) All 3-tuples of sym- 
bols selected repetitively from M, by requiring that each 
3-tuple is anti-symmetric with respect to its three entries, 
and in particular, if there are identical entries in a 3-tuple, 
setting the 3-tuple to be zero; such a 3-tuple is denoted by 

[ViVjVfe]. 

The two kinds of symbols must satisfy the following 
dimension-three constraints: 



IGP : For any five symbols Vi 1 , . . . , Vi 5 , 
[viiVijHv^v^Vig] - [v, 1 v, 3 ][v, ] v, i v, 5 ] 

+ [Vi 1 V j4 ][v i2 Vi 3 V i5 ] - [Vi 1 Vi 5 ]{v t2 Vi 3 V M } = 0. 



DB : For any six symbols 
[v 4l v i2 v 13 ][v i4 v 15 v i6 ] = 



[v^ViJ [v^Vig] [ViiVig] 
[Vi 2 Vi 4 ] [Vi 2 Vi 5 ] [Vi 2 Vi 6 ] 

[vi 3 Vi 4 ] [v i3 v i5 ] [v i3 v i6 ] 

(2.12) 



The inner-product bracket algebra is the commutative ring 
generated by the above two kinds of symbols, satisfying 
the symmetry requirements and the dimension-three con- 
straints. 

To include brackets of longer length, the concept quater- 
nionic bracket algebra or Clifford bracket algebra needs to 
be introduced. As explained before, there is no need to dis- 
tinguish between the two concepts in the setting of 3D or- 
thogonal geometry, so we simply call it bracket algebra. To 
distinguish from the concept of the same name arising from 
Grassmann-Cayley algebra |14] . we call that in [14] classical 
bracket algebra. 

Formally, besides the above 2-tuples and 3-tuples, a hier- 
archy of infinitely many new symbols can be defined: for 
any length I > 3, there are all Z-tuples of symbols selected 
repetitively from M, with the requirement that the first and 
the last equalities in Caianiello expansion (12. 9[) are satisfied; 
such an Z-tuple is denoted by [vj 1 Vj t • • • v,,]. By setting 
[1] — 1 (0-tuple) and [v^ — (1-tuple) for all i, we get a 
full hierarchy of new symbols marked by brackets, with ar- 
bitrary length I > 0. The new symbols together with their 
specific requirements, form a a commutative ring called the 
bracket algebra over 3D inner-product space. This is the the 
bottom-up approach to defining bracket algebra. The con- 
cept quaternionic product or Clifford product is not needed. 

3. VECTOR- VARIABLE POLYNOMIALS 
AND BRACKET POLYNOMIALS 

In this paper, we use (1) bold-faced digital numbers and 
bold-faced lower-case letters to denote vector variables, e.g., 
v, 1; (2) bold-faced upper-case letters to denote monic vector- 
variable monomials, e.g., A, X; (3) Roman-styled lower-case 
letters to denote polynomials, e.g., f,g; (4) Greek letters to 
denote K-coefficients, e.g., A, /x. 

Although the background is real orthogonal geometry, the 
algebraic manipulations under investigation are independent 
of the real field. In fact, only the following coefficients occur 
in computing: ±2 fe for k £ Z. We set the base field K to be 
of characteristic 7^ 2. 

Now start from quaternions. Let vi, V2, . . . , v n be symbols. 
What properties determine that the multilinear associative 
product among the symbols is the quaternionic one, and that 
these symbols represent vectors of a 3D real inner-product 
space with metric diag(— 1, — 1, — 1)? 7 gives a rather simple 
answer. 

The inner product of two vectors v;,Vj is a scalar, so it 
commutes with a third vector v^: [viVj]vfc = VfefviV,-]. For 
three vectors Vj, v,-, Vfc, since [viVjVfc] is a scalar, for a fourth 
vector V;, the commutativity [viVjVk]vi = v;[viVjVfc] holds. 
The two commutativities are all that characterize the equal- 



ity properties of the quaternionic product, besides the mul- 
tilinearity and associativity. 

Theorem 1. [7] Let vi, V2, . . . , v„ be n > 2 symbols. De- 
fine the product among them, denoted by juxtaposition of el- 
ements, as the WL-tensor product modulo the two-sided ideal 
generated by the following tensors: 

V2 



V3 
V4 



(Vj ® Vj + Vj ® Vi) ® v fc — Vfc IS 

(vj ® Vj ® v fe — Vfc ® Vj ® Vi) IS 

-V; ® (Vi ® V 3 ® Vfc - Vfc 



1 (Vj 1 
1 V; 



) Vj + Vj 
1 Vi), 



(3.1) 

for any i 7^ j 7^ k 7^ I in 1,2, ... ,n. Denote by Q or 
Q[[vi, . . . , v n ]] the WL-algebra defined by the above product 
and generated by the Vi. Denote by I[[vi, . . . , v n ]] the above 
ideal, and call it the syzygy ideal of Q. 



Denote 

K 2 := K({vi®Vi, Vj ® Vj + Vj ® v», \i ^ j}} 

»j Vfc — Vfc ® Vj 



K 3 := K 2 ({vi®Vj 



vi|i^j^fc}). 

(3.2) 



(1) Let l := [V1V2V3]. Then 1 / 0, and K 3 = K 2 (t)- 

(2) The Vi and iVj x Vj of the K-algebra Q span a 3D K2- 
space V 3 . For any k > 1, A(v J1 ■ ■ ■ v J2k+1 ) and iA(v jl ■ ■ ■ v J2k 
are both in V 3 . 

(3) The defined product is the Clifford product of the K2- 
Clifford algebra over V 3 . 

(4) 7/K = K3 = R, and the inner product of real space V 3 
induced from \2.1\l is definite, then the defined product is the 
quaternionic product. 



(5) Let Vfc = Vj 1 v ]2 ■ 
hold modulo I[[vi, V2, 



Then the following identities 



VitV fc ] = [Vfc] Vl , [ Vl V fe ] = [VfcVi] 



(3.3) 



The requirements in (4) distinguishing the quaternionic prod- 
uct from the Clifford product cannot be represented by equal- 
ities. So for symbolic manipulations of equalities, the quater- 
nionic product and the Clifford product cannot be distin- 
guished. The K-algebra Q is called the 3D vector-variable 
polynomial ring generated by vector variables Vi, and the 
product in it is called the vector-variable product. It is nei- 
ther the quaternionic product nor the Clifford product, but 
a more basic one. 

All the terminologies introduced earlier on quaternions and 
Clifford algebra are valid for Q. Besides, there are some 
additional terminologies for Q. Let vi -< V2 -< . . . -< v n 
be vector variables. A monic monomial of vector variables 
refers to the product of a repetitive permutation of some of 
the vector variables. For a monic monomial Vi 1 Vi 2 ■ ■ ■ Vi k , 
the leading variable refers to v^ , and the trailing variable 
refers to Vi k . The monomial is said to be non-descending if 
Vi 1 < Vi 2 ^ . . . X Vi k , and is said to be ascending if v^ -< 
Vi 2 -<...-< Vi k . The degree, or length, of the monomial is 
k. The lexicographic ordering among monomials is always 
assumed. 



A polynomial of vector variables is a K-linear combination 
of monic monomials. The leading term of a polynomial / is 
the term of highest order, denoted by lt(/). The degree of 
a polynomial is that of its leading term. The leading terms 
of all elements in a subset S of polynomials are denoted by 
It (S). When specifying the field K2 or K3, we can get the 
corresponding concepts quaternionic polynomial and Clifford 
polynomial. 

Fix a multiset of vector variables M composed of m > 3 
symbols vi, V2, • • • , v m . Let n be the number of different 
elements in M, where 3 < n < m. In the K-tensor al- 
gebra <2$(vi, V2, . . . , v n ) generated by the n symbols taken 
as vectors, a tensor monomial is up to coefficient the ten- 
sor product of finitely many such vectors. The K-tensor 
algebra over M, denoted by <^)[A4], is the K-subspace of 
<SK v i) V2 > • • • 1 v «) spanned by tensor monomials whose vec- 
tor variables by counting multiplicity are in M, equipped 
with the tensor product that is undefined if the result is no 
longer in ®[.M]. 

When the product among the elements in M is the vector- 
variable product, we have the corresponding concept of 3D 
vector-variable polynomial ring over multiset M, denoted 
by Q[yM]. Each element in Q[.M] is a polynomial whose 
multiset of vector variables in each term is a submultiset of 
M. The syzygy ideal X[M] of Q[M] is still defined by (f3~Tj) . 

The concepts of Grobner base and normal form are defined 
in Q[.M] just as in ®(vi, V2, . . . , v„) [9]. For two monomi- 
als hi,hi in vector variables, hi is said to be reduced with 
respect to /12, if hi is not a factor of hi, or hi is not a multi- 
plier of /12, i.e., there do not exist monomials I, r, including 
elements of K, such that hi = lh2r. For two polynomials 
/ and g, f is said to be reduced with respect to 5, if the 
leading term of / is reduced with respect to that of g. The 
term " non-reduced' means the opposite. 

Let /2, . . . , ft} be a set of vector- variable polynomi- 
als. Another set of vector-variable polynomials {gi, gi, . . ., 
gi} is said to be a reduced Grobner base of the ideal X := 
{h-,h, ■■■,fk) generated by the ft in Q[M], if (1) (31, • ■ • , 9l) = 
X, (2) the leading term of any element in I is a multiplier of 
the leading term of some gi, (3) the gi are pairwise reduced 
with respect to each other. 

The reduction of a polynomial / with respect to a reduced 
Grobner base <?i, <?2, • • • ,9l is the repetitive procedure of di- 
viding the highest-ordered non-reduced term L of / by some 
gi whose leading term is a factor of L, then updating / by 
replacing L with its remainder, until all terms of / are re- 
duced. The result is called the normal form of / with respect 
to the Grobner base. Two polynomials are equal if and only 
if they have identical normal forms. 

In [7J, two theorems are established for the Grobner base 
and normal forms of 3D vector-variable polynomials, one for 
the multilinear case where each element in multiset M has 
multiplicity 1, the other for the general case Q[[vi, . . . , v n ]]. 

Theorem 2. 7 Let Z[vi, . . . ,v n ] be the syzygy ideal of 
the multilinear polynomial ring Q[vi, . . . , v n ] in n different 
vector variables vi -< V2 -< . . . -< v„. 



(1) [Grobner base] The following are a reduced Grobner base 
ofX[vi, . . . , v n ]: for all 1 < ii < ii < . . . < ij < n, 

G3: [vi 3 Vi 2 Vj 1 ]-[vi I Vi 3 v< 3 ], and[\i 3 v il -v i2 ]-[vi 2 v i3 'Vi 1 ]; 

Gj: [vi 3 Vj 2 Vi 4 v i5 • ■ • v^.ViJ - [v i2 v l4 v l5 ■ • ■ v^-VijViJ, for 
all j > 3. 

(£) [Normal form] In a normal form, every term is up to 
coefficient of the form VYj v ai Vy 2 v Z2 • • • Vy fc v Z(c or ~V Yl V Z1 
■■■V Yk v Zk V Yk+1 , where 
(1) k>0, 

(u) v zl v Z2 • ■ • v Zfc is ascending, 

(Hi) every W ( is an ascending monomial of length > 1, 

(iv) VyjVy 2 ■ • ■ V Yk (or V Yl Vy 2 ■ • ■ Vy t+1 i/Vy t+1 occurs) 
is ascending, 

(v) for every i < k, ifv ti is the trailing variable of monomial 
V Yi , then v Zi -< v ti . 

Theorem 3. 7 LetX[[vi, . . . , v n ]] be the syzygy ideal of 
the polynomial ring Q[[vi, . . . , v„]] in n different vector vari- 
ables Vi -< V2 -< . . . -< v„ . 

(1) [Grobner base] The following are a reduced Grobner base 
ofX[[vi, . . . , v„]|: 

G3, Gj: for all 3 < j < 00, and all 1 < ii < 12 < is < ii < 
. . . < ij-i < ij < n; 

EG2: for all ii < ii, 

Vi 2 V 12 V n - Vi! Vjj, Vi 2 , Vi 2 Vi!Vi! -V^V^V^; 

EG j : [vi 3 Vj 2 Vi 3 Vi 4 • ■ • Vi.Vij] - [v i2 v i3 v l4 ■ ••v^v <1 v is ], for 

all 2 < j < 00, and all 1 < ii < «2 < *3 < *4 < < 

ij-i < ij < n. 

(2) [Normal form] In a normal form, every term is up to 
coefficient of the form V Yl v hl v zl V Y2 v h2 v Z2 ■ ■ -~V Yk v hk v Zk 
or V Yl v hl v zl ■ ■ ■ VY k v hh \ Zk V Yk+1 , where 

(I) k>0, 

(II) v zl v Z2 ■ ■ ■ v Zk is non-descending, 
(Hi) Vh 1 Vh 2 ' ' ' v h fc ts non-descending, 

(iv) every Vy ; is a non-descending monomial of length > 0, 

(v) V Yl v hl V Y2 v h2 ---V Yk v hk (orV Yl v hl ■ ■ ■ V Yk v hk V Yk+1 
i/Vy fc+1 occurs) is non-descending, 

(vi) for every i < k, y v Zi , 

(vii) for every i < k, if the length of~V Yi is nonzero, let wt i 
be the trailing variable of~V Yi , then vt 4 -< v^. 

For a general multiset M in which the different vector vari- 
ables are Vi, V2, . . . , v n , the Grobner base of the syzygy ideal 
Z[.M] is the restriction of the Grobner base of X[[vi, . . . , v n ]] 
to Q[jM], denoted by In Q[A4], a polynomial is said 

to be X-normal if its leading term is reduced with respect 
to the Grobner base. The procedure of deriving the normal 
form of a polynomial is called X-reduction. 

Now that any vector-variable polynomial has a normal form 
by X-reduction, so does a bracket polynomial when every 



bracket is expanded into two terms by definition. The result 
is complicated. 

Consider the following simple example: for a single bracket 

2[viV 2 ■ • ■ V m ] = V1V2 • ■ ■ V m + (-l)™V m • ■ • V2V1, 

the X-reduction goes as follows: for < j < m, if we define 
V m _j = Vj + iVj + 2 ■ • ■ Vm, then when m > 3, 

V ra + (-l) m V+, 

= ViV m _i + V m _ lVl - Vi(V m -i + (-ir^V^) 
= ViV m _i + V m _!V! - Vi(v 2 V m _ 2 + V m _ 2 V 2 ) 

+v lV2 (V m _ 2 + (-l) m - 2 Vl- 2 ) 

= V 2 V m _2Vi - ViV m _ 2 V 2 

+v lV2 (V m _ 2 + (-l) m - 2 V^_ 2 ). 

(3.4) 

From this recursive formula, we get for m > 1, 

or 1 1 1 + (zlT 
2[viv 2 ---v m J = viv 2 ---v m 

m — 1 

+ ^2 (-l) I+1 (vi • •• Vi- •■ V m )Vj. 

i=l 

(3.5) 

So the normal form of the simplest bracket [V1V2 • • • v m ] is 
composed of up to m terms. What is worse is that only 
when the terms are summed up can they represent a single 
algebraic invariant (the bracket), while missing a single term 
destroys the invariance of the whole expression. 

From the appearance, the bracket symbol only hides half of 
a binomial. There are much more behind this appearance. 
By definition, for a sequence A of a > 1 vector variables, 
its bracket is [A] = 2 _1 A + (-l)°2 _1 A t . Monomial A is 
called the representative of bracket [A]. Later on, when 
we write [A], we always assume that monomial A is the 
representative of the bracket. 

The definition of a bracket endows the symbol with the re- 
version symmetry (or equivalently, the conjugate symmetry) 
up to sign: [v^ . . . VjJ = (— l) a [v ia ■ ■ ■ ViJ. By (|3.3[) . the 
bracket symbol also has shift symmetry. So up to sign a 
bracket of a vector variables has the symmetry group D 2a 
(dihedral group). 

The above analysis is only for a single bracket. For bracket 
polynomials, there are a lot of polynomial identities, or syzy- 
gies, among them. These complexities justify the separation 
of bracket algebra from vector-variable polynomial ring in 
symbolic manipulations of algebraic invariants. Finding the 
Grobner base of the syzygies and then characterizing the 
normal forms of bracket polynomials are the main goal of 
this paper. 

The following are some terminology on brackets. The repre- 
sentative of a bracket polynomial is the the vector-variable 
polynomial whose terms are each the product of the coeffi- 
cient with the representatives of the bracket factors in the 
same term. The representative of a bracket polynomial is 
allowed to contain brackets. For example, [viV2[v3V4]] is 
taken as bracket binomial 2 _1 [viV2(v3V4 + V4V3)]; its repre- 
sentative is the content ViV2[v3V4] within the outer bracket. 



The lexicographic ordering of bracket polynomials is that 
of their representatives. The leading variable of a bracket 
refers to that of its representative. The leading term of a 
bracket polynomial is always under the lexicographic order- 
ing. For example, if Vi -< v 2 , then [viv 2 ] -< [v 2 Vi], and 
[viv 2 ][v 2 vi] -< [v 2 vi][viv 2 ]. 

The leader (or expanded leading term) of a bracket polyno- 
mial, refers to the leading term of the bracket polynomial 
when taken as a vector- variable one, i.e., the corresponding 
vector-variable polynomial obtained from expanding each 
bracket into two terms by definition. For example, the leader 
of bracket [A] refers to the one of higher order between 2 _1 A 
and (-l) a 2- 1 A t . 

Among all the bracket polynomials that are equal to the 
same bracket polynomial, there are two that have strong 
features: the first is the one whose representative is the low- 
est, the second is the one whose leader is the lowest. The 
second is unique but the first is not. To make the first unique 
we introduce the following concept. 

In Q[A4], where the number of elements in multiset Ml is m, 
a uni-bracket monomial refers to a single bracket of length 
m. A uni-bracket polynomial is a K-linear combination of 
uni-bracket monomials. All uni-bracket polynomials form a 
K-linear space, denoted by [Q][jM]. The K- linear space of 
the representatives of elements in [Q] [Mi] is just the K-linear 
space of degree-m vector-variable polynomials, denoted by 
Qm[jM]. Obviously, [Q][jM] is a linear subspace of Q m [.M]. 

When taken as a vector-variable polynomial, a uni-bracket 
is a binomial. In appearance, a uni-bracket is a monomial. 
To distinguish between the two understandings, we need a 
device to get rid of the bracket symbol and extract the rep- 
resentative of the uni-bracket. This can be done by taking 
[Q][A1] as the quotient of Q m [A1] modulo the ideal 

J[M]:=X[M] + [X][M], (3.6) 

where [X][A1] is composed of the vector parts of degree-m 
polynomials, i.e., the K-linear span of elements of the form 

R: v it v i2 . . .v im - (-l) m v lm ■ • ■ Vi 2 Vi x , (3.7) 

for all permutations of the m elements in Mi. 

The modulo- [X] [M] operation identifies a uni-bracket with 
its representative, or equivalently, identifies any degree-m 
vector-variable polynomial with the uni-bracket polynomial 
it serves as the representative. This operation retains the 
bracket symbols of all the brackets of length < m, while 
removing the bracket symbols from all brackets of length m. 

4. GROBNER BASE AND NORMAL FORM 
FOR MULTILINEAR UNI-BRACKET 
POLYNOMIALS 

From this section on, we use bold-faced digital numbers to 
denote vector variables, and use bold-faced capital letters to 
denote monic monomials of vector variables. 

In this section, the multiset Mi is composed of m > 3 differ- 
ent vector variables, and the modulo- [I] [M] operation is al- 
ways assumed. Then [Q][jM] and Q m [jM] are identical, and 



a uni-bracket no longer has the outer bracket symbol. Ideal 
[X][.M] is called the uni-bracket removal ideal in Q[A4], and 
ideal J[M] is called the syzygy ideal of [Q][A4] in (g)[yW]. 
Below we compute the Grobner base of J\M\ and charac- 
terize the normal forms of uni-bracket polynomials. 

The following are elements of ^[TVf]: for all monomials A 
to F such that A1B, 1C, 1D2, 1E2F are of length m and 
E has length e > 0, 

Si: A1B-1BA, 

Sf : the Si in which A1B is X-normal, 

Ri : 2(1[C]) = 1C - (-l) m lC + , (4.1) 

R12 : 2(1[D2]) = 1D2- (-l) m 12D + , 

R12M : 2(1[E2[F]]) = 1E2[F] - (-l) e 12E t [F]. 



Since the reduced Grobner base of X[M] is 5[A4], we only 
need to consider the elements of type R in i7[.M], as they 
span [X][M]. By 

A1B - (-l) m B + lA t 
= (AlB-lBA) + (lBA-(-l) m lA t B t ) (4.2) 
-(-l) m (B t lA t - lA+Bt), 

we get 



Lemma 4. R is a subset of the ideal (Si,Ri). For any 
type-R element f but not of type Ri, lt(/) G lt(Si) . 



Lemma 5. Si is a subset of (Sf ,Ri) +1[M]. For any 
type-S! element f but not of type Sf , lt(/) G \t{I[M\). 



Proof. Consider a general type-Si element / — A1B — 
1BA, where A is not empty. If A1B is X-normal, then / G 
Sf. If X-reductions are carried out to A, B, say A = A N +X 
and B = B^ + X, where A^B^ are both X-normal, then 
/ = A^IB^-IBV+I, and A^IB^ is the new leading 
term. So by X-reductions we can assume that both A and 
B are X-normal. 

Assume that A1B is not X-normal. Further assume that 
any type- Si element g -< / is in (Sf , Ri) + X[.M]. We prove 
the conclusion for / by reduction on the order. There are 
three possibilities to apply Grobner base elements of X[A4] 
to make reduction to A1B: (i) G3 at the end of Al, (ii) Gi 
for i > 3 at the end of Al, (iii) G3 on 1 and two variables 
from A, B respectively. 

Case (i). Let A1B = CuvlB, where u >- v >- 1. Then 
/ = CvulB + Cl(uv + vu)B lBCuv 

induction _ 1BC(vu+uv) + l(uv + vu)BC 
I 0. 

Case (ii). Let A1B = CuvDIB, where 1 -< v -< u -< all 
variables of D. Let the length of D be d > 0. Let the lengths 



of B, C be b, c respectively. Then b + c = m — d — 3. 

/ = C(vDl + (-l) d lD t v)uB-(-l) d CulD t vB 

-lBCuvD 

ind = tion l{ (uBC - BCu)vD + (-l) d Dt v (uBC - BCu)} 

= x (-l) d lD t v{(-l) m - d (C t B t u - uCtBt) 

+(uBC-BCu)} 
= (-l) d lD t v{-((-l) 6+c C t B t + BC)u 

+u((-l) f,+c C t B t + BC)} 

= 0. 

Case (iii). Let A1B = CulvD, where u >- v >- 1. Let the 
lengths of C, D be c, d respectively. Then c + d — m — 3. 

/ = -CluvD + Cv(ul + lu)D-lvDCu 
induction 1 (_ m2 c + DCvu + uDCv - vDCu) 
l{-((-l) m C t D t -DC)vu 

+v((-l) m C t D t -DC)u} 

I 0. 

□ 



Lemma 6. Ri is a subset of the ideal (R12) +X[A^]. For 
any type-Ri element f but not of type R12, lt(/) G lt(X[7V(]). 

Proof. Let / = 2(1 [E2F]) be a general clement of type 
Ri. Then / = 2(1[FE2]), and the latter is in R i2 . □ 

Theorem 7. Let M = {1, 2, . . . , m} be m > 2 different 
symbols, where 1 -< 2 -<...-< m, and let Q[M] be the 
vector-variable polynomial ring over M. Let [Q][M] be the 
space of uni-bracket polynomials, and let J[M\ be its syzygy 
ideal in ®[M]. 

(1) [Grobner base] The following are a reduced Grobner base 
of J[M]: 

G[M\: Gifor all3<i< m; 

Sf : A1B — IB A, where A is an ascending sequence of 
length a > 0, B is of length in — a — 1 > 0, and A1B 
ts X-normal; 

Rf 2 [j]: for all 1 < j < m=±, l[A2][Y 2 z 2 ][Y 3 z 3 ] • • • [Y jZj ], 
where (i) A is an ascending sequence of length a > 0; 

(ii) when j = 1, then Rf 2 [l] = 1[A2]; 

(iii) each Yj is a non-empty ascending sequence, and 
each Zj is a variable such that ZiYi is ascending; 

(iv) IA2Y2Z2 • • • YjZj is X-normal and length-m. 

(2) [Normal form] In a normal form, every term is X-normal, 
and is up to coefficient of one of the following forms: 

(I) 12C, where C is X-normal; 

(II) IA2Y2Z2Y3Z3 ■ • • Yfc_iZfc_iYfc, where k > 2, A and 
the Yj are each a non-empty ascending sequence, and each 
Zj -< ti, the latter being the trailing variable o/Yj ; - 



(Ill) IA2Y2Z2Y3Z3 ■ ■ ■ Y fc z fc , where k > 2, A and the Yi, Zj 
are as in (77j, ana" for some 2 < i < k, if 1, is tte leading 
variable ofYi, then L, -< Zj. 

Proof. There are several steps. 

Step i. We need to prove that R12 is a subset of the ideal 
(R12} +I[-M]. Once this is done, then since the leader of any 
element of Ky2 is I-normal and cannot be cancelled by the 
leader of any other element of RfJ, , the Gi and [j] form 
a reduced Grobner base of (R12} + I[.M]. By this and the 
previous three lemmas, we get conclusion (1) of the theorem. 

Once conclusion (1) holds, then any I-normal monomial of 
length m with leading variable 1 is the representative of a 
uni-braeket in normal form if and only if it is not the leader 
of an Rf^-typed element. Conclusion (2) follows. 

Step 2. The idea of proving the statement in Step 1 is to use 
C?[.M] to decrease the order of the leader of every element 
of type R12 or Ri2[*], at the same time keep the reduction 
result to be within the K-linear space spanned by elements 
of type R12 or Ri2[*]. Then ultimately all the leaders of 
these elements become I-normal. 

We start with the I-reduction on the leading term of a gen- 
eral Ri2-typed element / = 1A2 - (-l) m 12A + , where the 
length of A is m — 2. If by I-reduction, A = A^ + I, then 
/ = 1A N 2 - (-l) m 12A Jvt +1. So we can assume that A 
is I-normal. 

If A2 is I-normal, then / is just R^[l]. When A2 is not 
I-normal, if A2 is non-reduced with respect to G3, let A2 = 
Buv2, where u >- v >- 2, then 

1Buv2 = -1Bvu2 + 2{lB2[uv]}. (4.3) 

The result consists of the leading terms of one Ri2-typed 
element and one Ri2[*]-typed element. The leaders of both 
terms are lower than /. 

If A2 is non-reduced with respect to Gi for some i > 3, let 
A2 = BuvC2, where u >- v >- 2, and uC is ascending, and 
the length of C is c > 0. Then 

1BuvC2 = lB(vC2 + (-l) c 2C t v)u- (-l) c lBu2C t v 
= lByC2u + lB2((-l) c C t vu + uvC) 
!B2uvC + lBu2(-(-l) c C t v + vC) 
!Bu2vC 

= 2{lB2[uvC]} + 2{lBu2[vC]} 1BvCu2, 

(4.4) 

The result consists of the leading terms of two Ri2-typed 
elements and one Ri2[*]-typed element. The leaders of the 
three terms are lower than /. 

By ()4.3[) an d (|4.4[) . a monomial that is non-reduced with 
respect to a Gi for some i > 3 must contain a subsequence 
of the form uDv, where (i) the length of D is d > 0; (ii) 
u y v; (iii) if Id is the leading variable of D, then u >- Id; 
(iv) if D contains more than one variable, then Id >- v. (14.31) 
and (14.41) can be written in the following unified form: 

uDv G( = 2) 2(uv[D]) + 2(v[uD]) - Duv. (4.5) 



It is called the fundamental I-reduction formula. 

Step 3. Consider I-reductions on the leader of a general 
Ri 2 [*]-typed element / = 1A2[B] - (-l) a 12A f [B], where 
the length of B is b > 0. 

Since [B] = (— 1) [B'], henceforth we assume that in any 
Ri2[*]-typed element to be normalized, the leading variable 
of any bracket has higher order than the trailing variable 
of the bracket. Then the leader of the bracket is always its 
representative. 

In this step, we consider I-reduction to the representative 
B of [B]. Let [B] = [CuDvE]. Substituting (g3| into it, 
we get 

[CuDvE] = 2[CuvE][D] +2[CvE][uD] - [CDuvE]. (4.6) 

Step 4- Consider I-reductions on the leader 1A2B of / = 
1A2[B] - (-l) a 12A t [B] involving both the tail part of A 
and the head part of B, where the leading variable of B is 
assumed to be higher than the trailing variable. 

As 2 is lower than any element of A,B, the only possible 
reduction is by G3. Let 1A2B = lCa2bD, where C may 
be empty but D is not. Assume a y b >- to, where to 
is the trailing variable of D. Let the length of D be d. It 
is easy to prove that applying G3 to a2b in vector-variable 
binomial lCa2[bD] is equivalent to the following absorption 
of bracket: 

lCa2[bD] = lC[bD]a2 

= 2~ 1 (lCbDa2) + (-l) d 2- 1 (lCD t ba2). 

(4.7) 

Each term in the result is a leading term of an Ri2-typed 
element lower than /. 

Step 5. In Step 3, we have seen that a single bracket after 
I-reduction, may be split into two brackets. The split can 
continue and we gradually get expressions of the form 

R 12 [j] : 1E2[F 2 ][F 3 ] • • ■ [F f ] - (-lJ'UE* [F 2 ] [F 3 ] ■ ■ ■ [P f ], 

(4.8) 

where the length of E is e > 0, and the length of 1E2F 2 F 3 
• ■ • Fj is m. R12L7] is a K-linear combination of elements of 
type R12 [*] if all but one bracket are each expanded into two 
terms. 

Consider I-reductions of R12L7] involving more than two 
bracket factors, and I-reductions involving 1E2 and more 
than one bracket factor. Since 2 is lower than all elements of 
E and the Fj, G3 is the only possible Grobner base element 
that may apply to 2 and its neighbors on both sides simul- 
taneously. G3 can involve only [F2] among the brackets. 

In [F2][F 3 ] • ■ • [Fj], only Gi where i > 3 can involve more 
than two brackets. However, since Gi is of the form uDv 
where D is ascending, if the product of the leaders of three 
brackets is non-reduced with respect to some Gi, then the 
middle bracket must be composed of a subsequence of D 
of length > 2, contradicting with the assumption that the 
leading variable in the middle bracket be higher than the 
trailing variable. 



So each Z-reduction of R12L7] by a single Gi where i > 3, 
can involve at most two bracket factors, or the 1E2 and one 
bracket factor. 

Step 6. Consider X-reductions on [Fi] [F2] , where the leading 
variable in each bracket is higher than the trailing variable. 
If the leading variable If! of Fi is higher than the leading 
variable lp 2 of F2, then an X-reduction commuting the two 
brackets reduces the order of their product. Below we always 
assume \? 1 -< 1f 2 - 

For G3, there are two possibilities to involve both Fi,F 2 
in the leader F1F2: two variables at the end of Fi and the 
third at the beginning of F2, or one variable at the end of 
Fi and the other two at the beginning of F2. The latter is 
impossible because 1f x -< If 2 - For Gi where i > 3, there 
are also two possibilities: two variables at the end of Fi and 
the rest at the beginning of F2, or one variable at the end 
of Fi and the rest at the beginning of Fa. The latter is also 
impossible due to If 1 -< 1f 2 - 

Case G3. Let [Fi][F 2 ] = [Buv][wCd], where u >- v, and 
u >- w >- d. Let the length of C be c > 0, and let the 
leading variable of Bu be 1. Then w y 1 y v. Applying 
G3 to uvw is equivalent to the following absorption of the 
second bracket: 

[Buv][wCd] 

= [B[wCd]uv] (4.9) 
= 2~ 1 [BwCduv] + (-l) c 2^ 1 [BdC t wuv]. 

The leader of each bracket monomial in the result has lower 
order than the leader of [Fi][F 2 ]. 

Case Gi. Let [Fi][F 2 ] = [Buv][aDwC], where aD is as- 
cending, and a >- u >- v >- w. Let the lengths of B, C, D be 
b, c, d respectively. Let the leading variable of Bu be 1, and 
let the trailing variable of wC be t. Then a >- 1 >- v and 
a y t. Applying G(d + 4) to uvaDw, we get 

4 [Buv] [aDwC] 

= Bu vaDw C + (-l) c+d B uyCjjw D t a 

+ (-l) b vuB 1 aDwC + (-l) 6+c+[i vuB 1 'C t wD t a 

= (vaDw - (-l) d wD + av)BuC + (-l^BuwD+avC 
+(-l) d ((-l) c vC t w + wCv)BuD t a 
-( l^BuwCvD+a 

+ (-l) 6 v(aDwC + (-l) c+d C t wD t a)uB t 

= v{aDw(BuC - (-l) i ' +c C t uB t ) 
+ (-l)'aDw(C + (-lyC^uB^ 
+ (-l) c+d C t w(BuD t a+ (-lJ^aDuB') 
- (-l^C^aD - (-l) d D t a)uBt} 
+(-l) d w{-D t av(BuC - (-l) i,+e C t uB t ) 

- (-l^+TjWctuB 1 " - (-l) f,+d CvaDuB t 
+ Cv(BuD t a + f-^+'aDuB*)} 
+ (-l) d Buw(D t avC - CvDta) 

I 4([vaDw][BuC] + (-l) d [wCv] [BuDta]) 
+(-l) ft v(aDC + (-l) c+d D t aC t )wuB t 
-(-l) b w(CvaD + (-l) c+d D t avC t )uB+ 
+ (-l) d Buw(D t avC - CvD+a) 



= 4([vaDw][BuC] + (-l) d [wCv][BuD t a]) 
+(-l) ft 2{v(aD[C] - (-l)°[aD]C t )wuB t } 
-4 [Buw][CvaD] + Buw{C(aD - (-l) d D t a) 

+ (-l) d D t a(C + (-l) c C t )}v 

= 4{[vaDw][BuC] + (-l) d [wCv] [BuD f a] 

- [Buw][CvaD] + [Buw(C[aD] + (-l) d D t a[C])v]} 

= 4{[vaDw][BuC] + (-l) d [wCv] [BuD'a] 

- [Buw][CvaD] + [BuwCv][aD] 
+ (-l) d [BuwD t av][C]}. 

(4.10) 

The leader of each bracket monomial in the result has lower 
order than the leader of [Fi][F 2 ]. 

Step 7. Consider a bracket of the form h = [aiBiCia 2 B2C2 
•••afcBfcCfe], where (1) k > 1, (2) a^ y a for every i, (3) 
ai >- Cj for all 1 < j < k. Let the length of Bi be bi. 

When k = 2, 

[aiBiCia 2 B 2 c 2 ] - 2[aiBiCi][a 2 B 2 C2] 
= 2- 1 {(-l) bl+i,2 c 2 Bta2CiB t 1 a 1 - (-l) &2 aiBiCiC2B|a2 
-(-l) 6l ciB£aia3B a ca - (-l) 6l+b2 ciBjaiC2Bja2} 

= -(-lJ^lcaBjaaaiBiCi], 
The leader in the result has lower order than h. 

For k > 2, 

[aiBiCi ■ • ■ a fe B fc c fe ] 

-2[aiBiCi ■ ■ • a fc _iB fc _iC fc _i][a fc B fc c fe ] 

= -(-l) 6fc [c fc B^.a fc (aiBiCi • ■ • afc_iB fc _iCfe_i)], 

and for [aiBiCi • • • afc-iB^-iCfe-i], the split into [aiBiCi 
■ ■ ■ afc_ 2 Bfc_ 2 Cfc_2][afc_iBfc_iCfc_i] can continue. In the end, 
we get 

[aiBiCi ■ ■ ■ a fc B fe c fc ] = 2 fc_1 [aiBiCi] • ■ ■ [a fc B fc c fe ] + g, 

(4.11) 

where g is a bracket polynomial whose leader has lower order 
than h. 

(14. lip can be used to split a long bracket whose representa- 
tive is X-normal. It can also be used in the converse direc- 
tion, to concatenate short brackets into a long one. 

Step 8. So far we have proved that for any Ri2-typed element 
1[E2] or any Ri 2 [j]-typed element 1[A2][F 2 ][F 3 ] ■ ■ • [Fj], as 
long as the leader is not I-normal, I-reductions can always 
be carried out to change 1E2 or lA2[F2][Fs] ■ ■ • [Fj] into 
the following form: 

1T = J2 \ a lE a 2 + wlA^tPftlpr^] ■ ■ ■ [F, 8j ] 

a 

+ ^> 7 12[D 71 ][D 72 ]---[D 7 J, 

(4.12) 

where the leading variable in each bracket is higher than the 
trailing variable, the 1E Q 2 and lAg2Fa 2 F / 3 3 • • • F^ are all 
I-normal. 

Since any I-normal form is of type Y1Z1 . . . Y^Zk or Y1Z1 
. . . YfcZfcYfc+i, it must be that 



(i) lE a 2 = Yizi, i.e., E Q = 34 ■ ■ ■ m. 



■ F 0j = Yizi 


...Y 3 z 


1A 8 2 = 


Yi Zl , 


Ffe = 


Y 2 z 2 , 







(i) is obvious. In (ii), the trailing variable of each bracket 
must be some z;. If an is Y^z^ . . . Y^+pZ^+p for some 
p > 0, for h < s < h + p, let the leading variable of Y s be 
la, then \ h y z h+p y z h+p -i y . . . y z h . By 1)4. lip . [F^J is 
split into 2 p \YhZh] ■ ■ ■ \Y h+ P ^h+ P ] plus some bracket mono- 
mials of lower leader. Then I-reductions continue to the 
terms involving such bracket monomials. Ultimately each 
bracket is of the form [YjZj]. 

In (|433, lE a 2 = (-l) m l(E a 2) t by Rf 2 [l], 1A S 2 [Fp 2 ] • • • 
[E>.] = (-l) m l(A^2[F fe ] ■ • ■ [F, 8j ])t by Rf 2 [j], and 12[D 71 ] 
■ ■ ■ [D 7fc ] = (-1)-1(2[D T1 ] ■ • ■ [D 7 J)t by X[M]. So IT - 

(-l) m lT t is reduced to zero by X[M] and Rf 2 . □ 

5. SQUARE-FREE VECTOR- VARIABLE 
POLYNOMIAL RING 

When A4 is a general multiset of vector variables, a square 
in Q[jM] refers to the product of a vector with itself. Denote 
vf := VjVi. It commutes with everything in Q[.M]. 

Proposition 8. In Q[M], let f = gvi h be a multiplier of 
vf. If f is not I -normal, then by doing I -reduction to g,h, 
together with rearranging the position of~v\ in each term, f 
can become X-normal. 

Proof. Suppose g, h are X-normal. There are three cases 
for f to be non-reduced with respect to the Grobner base 
G[M\. 

(1) If / contains as a factor the leader of Gk for k > 3, or 
EGj for j > 1 involving both of vf , then vf is preserved by 
the reduction with the Grobner base element. 

(2) If gvi is non-reduced, then switch the element of <5[.M] 
with respect to which gvi is non-reduced: 

Case EG2: Let gvi — AuViVi or Auuvj, where u >- V;. 
Then Auv^vf h — AvfuVi/i or Auirvf h = Avf uu/i. 
Case G3: Let gvi — Auwv; where u >- w and u >- Vj. Then 
Auwvf h = Avfuw/i. 

Case Gfc for k > 3 or EGj for j > 2: Let gvi = CuwDv^ 
where u >- w >- Vj. Then CuwDvf/i = CvfuwD/i. 

(3) If Vih is non-reduced, then switch the element of <5[.M] 
with respect to which v^/i is non-reduced: 

Case EG2: Let v^/i = v^v^zB or v^zzB, where Vj y z. 
Then gvfv^zB = gv^zvf B or gvfzzB = grzzvfB. 
Case G3: Let Vih = v^yzB where Vj y y and Vj >- z. Then 
flvfyzB = ffyzvfB. 

Case Gfc for fc > 3 or EGj for j > 2: Let Vih = v t yCzD 
where Vi y y >- z. Then pvfyCzD = gyvf CzD. 



In all the cases, the order of / is decreased while preserving 
V; . By induction on the order we get the conclusion. 

Proposition|H]suggests a "square-free normalization" of vector- 
variable polynomials, by moving all squares to a set free of 
any reduction operation, and maintaining the set of squares 
in a normal form. 

In a vector- variable monomial, let the set of squares be sep- 
arated from the remainder of the monomial by a symbol 
such that all elements on the right side of the symbol 
are squares. Two things need to be established before such 
a symbol can be used in algebraic manipulations: (1) alge- 
braic structure of the new symbolic system, (2) connection 
with the canonical system based on V2, V3, V4. 

Let 5 be a commutative monoid. All elements in S span a K- 
vector space whose dimension equals the number of elements 
in 5. The product in the vector space is the multilinear 
extension of the product in 5. The vector space equipped 
with this product forms a commutative K-algebra, called the 
K-algebra extension of monoid S, denoted by K<S. 

For a K-algebra A, when 5 is a subset of the center of A, 
then A is not only a module over the ring KS, but a multi- 
linear algebra over K<S, called a M.S -algebra. 

For K-tensor algebra (^)[A4], let 

® D {M] :=®[M]/(V2). (5.1) 

It is easy to see that when setting S to be generated by 
elements of the form vf := Vi ® Vi, for all Vi 6 M, then 
<^) D [jM] is a IKS-algebra, called the KS -tensor algebra over 
multiset M, or the square-free tensor algebra over M. The 
product in (^^Al] is induced from the tensor product. For 
brevity we still denote the product by but denote the 
commutative product in 5 by juxtaposition of elements. 

In D [A4], for all q € 0|yW] and s G 5, we introduce the 
notations 

qDs := q<S>s e®[M], 

qU := q &®[M), (5.2) 

□s := s G S. 

Then 

(giDsi) ® (g 2 Ds 2 ) := qi ® g 2 Dsis 2 . (5.3) 

Formally, an element q £ Q[A4] is taken as (jOl, and an 
element s £ IKS is taken as IDs. In other words, factor Dl 
(or □) in gDl (or q\3) is usually omitted. So 

qDs = q ® (Da) = (Da) ® q, (5.4) 

and 

Dst = (Ds) ® (Dt) = (Dt) ® (Us). (5.5) 

That S is generated by squares can be succinctly expressed 
by the following identity: 

Vi ® Vi = Ovf. (5.6) 

making left multiplication with / and right multiplication 
with g on both sides of the identity, we get f®(vi®Vi)®g = 



/ Cg> <pvf . It includes V2 as a special case. 

The degree, or length, of a monomial in <2> n [.M] is the degree 
of the monomial when taken as an element in (^)[.M]. The 
left degree or left length of a monomial refers to the degree 
of the monomial on the left side of the square symbol. For 
a monomial / £ <^) D [A1], its canonical form in <^)[A4] is 
denned to be the monomial of lowest lexicographic order 
among all monomials equal to / modulo V2. The order of 
/ is that of its canonical form. This ordering is still called 
the lexicographic ordering. 

The canonical form of f = v^, <g> Vi, ® • • • ® v., □v? ri v? r2 

* 1 z K Jl 72 

• ■ ■ v J( r ' , where -< v J2 -<.-.-< Vj, , can be obtained as 
follows: 

1. Set g = v» t ® Vi 2 ® • • • ® Vi,.. 

2. For p from 1 to Z, let v; t be the first variable in the 
sequence of g such that v^ t >- v Jp X v^^. Insert 

Vj p ® Vj p g) • • • <g) Vj p 
2r p 

to the position before v^ t in g, and update g. 

3. Output g. 

The vector-variable polynomial ring Q[.M] when taken as 
the quotient of (^) [Ai] modulo the two-sided ideal Z D [A4] 
generated by V3, V4, is a KS-algebra, called the square- 
free polynomial ring, denoted by Q D [A1]. The product in 
Q D [A1] is still denoted by juxtaposition of elements. I D [A4] 
is called the syzygy ideal of Q D [.M]. 

Theorem [3] has the following square- free version for multiset 
M: 

Theorem 9. Let Ai be a multiset of m symbols, among 
which n are different ones: vi -< V2 -< . . . -< v n , and let 
X n [Ai] be the syzygy ideal of the square-free polynomial ring 
Q a [M]. 

(1) [Grobner base] The following are a reduced Grobner base 
ofX n [M]: G3, and 

Gj : for all 3 < j < n + 1, and i\ < ii < . . . < ij, 

[Vi 3 Vi 2 Vi 4 Vj 5 ■ • ■ VijVjJ - [v, 2 V i4 Vj 5 ■ ■ • Vij-V^Vig], 

EGA;/ for all 3 < k < n + 1, and ii < 12 < . . . < %k, 

[v i3 v i2 v 13 v l4 ■ ■ • Vi fc ViJ - [v 12 v i3 v 14 ■ ■ • v ifc v n v !3 ]. 

The above Grobner base is denoted by Q U [M]. 

(2) [Normal form] In a normal form, every term is up to 
coefficient of the form Vyi v zl Vy 2 v Z2 • ■ ■ V^v^ds or 
V yi v zi • ■ • Vy,v Zt Vy HI Ds, where 

(i) k > 0, 

(ii) v zl v Z2 • ■ • v Zfc ts non-descending, 

(Hi) every is an ascending monomial of length > 0, 

(iv) V Yi Vy 2 ■ ■■ Vy k (orV Yl --- Vy t Vy i+1 i/Vr H1 occurs) 
is non-descending, 

(v) for every i < k, let v ti be the trailing variable of Vy ; , 



then vt; y v Zi , 

(vi) s is either 1 or the product of several squares. 

In Q n [A4] , a polynomial is said to be TP -normal if its leading 
term is reduced with respect to the Grobner base C/ D [.M]. 

A monic square-free uni-bracket monomial is of the form 
[A]Ds, where A is either 1 or a monomial of length > 1, 
s is a product of squares, and the length of As is m. A 
square-free uni-bracket polynomial is a K-linear combination 
of square-free uni-bracket monomial. The space of square- 
free uni-bracket polynomials is denoted by [Q] D [A1]. 

The K-linear space of degree-m square-free polynomials is 
denoted by Qml-M]. The space [Q] D [A4] can be taken as 
the quotient of Q^JAI] modulo the ideal 

J a [M]~I a [M] + [l] n [M], (5.7) 

where [I] n [7V(] is composed of the vector parts of degree-m 
square- free polynomials, i.e., the K-linear span of elements 
of the form 

R D : AOs - (-l) a A t Ds, (5.8) 

where A is a monomial of length a > and contains no 
square, s is a product of squares, and the length of As is m. 

The modulo- [X] n [Ai] operation identifies a square-free uni- 
bracket with its representative. It removes the outer bracket 
symbol on the left side of the "□" symbol from every square- 
free uni-bracket, disregarding the length of the bracket. Ideal 
[X] D [M] is called the uni-bracket removal ideal in Q^fAl], 
and ideal is called the syzygy ideal of [Q] D [A4] in [Ai]. 

6. GROBNER BASE AND NORMAL FORM 
FOR UNI-BRACKET POLYNOMIALS 

In this section, we extend Theorem [7] to the case of general 
multiset Ai with m > 3 different vector variables. The 
modulo- [X] ]Ai] operation is always assumed, i.e., ]Q!f i ]Ai] 
and Qm[A-1] are identical, and a uni-bracket does not have 
the outer bracket symbol on the left side of the "□" symbol. 

The following are elements of .^[A^]: 



R D (fc) : 


(K-(-l) fc Kt)Ds, 


Si(fc) : 


(iAbiB - biBiAps, for i / bi 


Sio(fc) : 


(iCbi - biiCps, fori^bi, 


Sn(fc) : 


(biAbiB - biBbiAps, 


S qi (fc) : 


(biCbi - CDbfps, 


SqiW(fc) : 


(biEbi - EDbf)[F]Ds, 


Ri(fc) : 


bi[D]Ds, 


Rn(fc) : 


bi[Cbi]Ds, 


Ril [*](*) : 


bi[Ebi][F]Ds, 


Rw(fc) : 


bi[Cb 2 ]Ds, 


R12 [*](*:) : 


bi[Eb 2 ][F]Ds, 



where 

(a) s has length m~k, and the left length of each expression 
is k > 0; 

(b) bi , i>2 are respectively the variables of the lowest order 
and the second lowest order on the left side of the square 
symbol; 



(c) in Si, either A or B can be empty, while in Sn, both A 
and B are non-empty; 

(d) in Ri, D is either empty (i.e., D = 1), or of length > 1; 

(e) in Sio and Rn, C is non-empty; 

(f) in Sq 1 [*](fc) and Rn[*], E, F are non-empty, and F does 
not contain bi; 

(g) in R12, C is non-empty and does not contain bi; 

(h) in R12 [*] , E, F are non-empty and do not contain bi , 
and F does not contain b2. 

In bracket [1A1], we have [1A1] = [Apl 2 . That the lead- 
ing variable has higher order than the trailing variable is 
always possible. This is taken as a postulate for all the 
brackets in (|6.1|l . 

Consider a general element / = (K — (— l) fc K t )Ds of type 
R D (fc): 

1. If bi occurs in K both as the leading variable and trailing 
variable, then / G (Sq 1 (fc),R D (fc - 2)). 

2. If bi occurs in K at only one end, then / G (Ri(fc), Si(fc)}. 

3. If bi occurs at the interior of K, set K = AbiB, where 
A, B are both non-empty, and bi does not occur at any end 
of A or B. By g^J), / G (Ri(fc), Si(fc)). 

By induction on k, we get 



Lemma 10. 

R°(fc) C + (SQaW) + <Ri(A)»- 

h<k 



g n [M]: Gi for all3<i< n; EGj for all 3 < j < n + 1. 

S?: (AbiB — biBA)Ds, where the length of each term is 
m, bi is the variable of the lowest order on the left 
side of the square symbol, A is a non-empty ascending 
sequence not containing hi, and AbiB is X n -normal; 

bi[Yibi][Y 2 bi] • ■ • [Y J bi][Y J+ iz J+ i][Y J+2 z J+2 ] 

■ ■ ■ [Yi+iz j+; ps, 

(6.3) 

where (i) j, I > 0, and the length of each term is m, 

(u) bi is the variable of the lowest order on the left side of 

the square symbol, 

(Hi) each Yi is a non-empty ascending sequence not con- 
taining bi, 

(iv) for all j + 1 < i < j + 1, z; 7^ bi, and ZiY % is ascending; 

(v) biYibi • • • YjbiYj+iZj+i • • • Yj+iZj+i is X D -normal. 

(2) [Normal form] In a normal form, every term is X n - 
normal, and is up to coefficient of one of the following forms, 
where statements (ii), (iii) on bi and the Yi are still valid: 

(I) biYibi ■ • ■ YjbiYj+iZj+i ■ • ■ Y j+ iz j+ iY j+ i +1 ns, where 
j, I > 0, and for all j + 1 < i < j + 1, z; -< ti, the latter being 
the trailing variable ofYi; 

(II) biYibi ■ ■ ■ YjbiY J+ iz J+ i • ■ • Y J+i z j+; Ds, where j,k > 
but j + k > 0, each Zi -< ti, but for some 1 < h < I, if\j+h 
is the leading variable ofYj+h, then lj+h ^ z j+h- 



In Ri(fc), when D contains bi, let D = AbiB, where the 
lengths of A, B are respectively a, k — a — 2, then 

bi[AbiBps = bi[BAbi]Ds 

= (biBAbi - (-l) fc A t B t )Db?s 

= (biBAbi - BADbfps (6.2) 

+(BA - (-l) fc A t B t )DbiS 
G <S qi (fc),R n (fc-2)>. 

By induction on k, we get that both bi[AbiB]Ds and 
bi[BAbi]Ds G Rn(fc) are equivalent to Sq x (fc): (biBAbi — 
BADbi)Ds in the sense that their difference is in the ideal 
E h < fc - 2 «Si (h)) + (Ru(ft)> + <R U (A)» + <Ri(l)) +T°[M]. 



Lemma 11. R D (fc) is a subset of the ideal 

X)«Si(A)> + <Rix(ft)> + (R12W)) + <Ri(l)) +T D [M]. 

h<k 

Theorem 12. Let M be a multiset of m > 2 symbols, 
among which n > 2 are different ones: 1 -< 2 -<...-< n, 
and let Q n [M] be the square-free vector-variable polynomial 
ring over Ai. Let [Q] n [A4] be the space of square-free uni- 
bracket polynomials, and let J~ D [A4] be its syzygy ideal in 
® a [M}. 

(1) [Grobner base] The following are a reduced Grobner base 
ofJ a [M], denoted by BQ[M\: 



Remark. The set of Ri can be replaced by the following 
three sets of degree-m polynomials: for all j > 0, I > 0, 

R?[0,0] : bids, 

Sq°[j,l] : (biYibi - YiDb?)[Y 2 bi] ■ • ■ [Y,bi][Y J+ iz,+i] 

• • • [Yj+(Zj+ips, 

B&V, I] : bi[Yib2][Y 2 b 2 ] ■ • ■ [Y J b2][Y J+ iz J+ i][Y J+2 z J+ 2] 

■ ■ ■ [Yj+tzj+ips. 

(6.4) 

In R?2 [j, I] i Y j+i, . . . , Yj + i do not contain b2, and the Zi >- 
b2. The replacement has no effect upon the normal forms. 

Proof. There are several steps. 

Step 1. We need to prove by induction on k that Si(fc), 
Rn(fc), Ri2(fc) are all in the ideal 

Y, <S?,R?M)+I D [M], (6.5) 

left length <fc 

where the asterisk stands for the (j, I). 

Once this is done, then since the leader of any element of 
type S? or R?[*] is I D -normal and cannot be cancelled by 
the leader of any other element of type S? or R? [*], the 
Si,Ri [*] and (J D [A1] must be a reduced Grobner base of 
(Si,Ru,Ri2,Ri(l)) +l n [M}. By Lemma [Til this ideal is 
just J a [M]. This proves conclusion (1), and conclusion (2) 
follows. 



Step 2. Si(fe) requires k > 1; Rn(fc) and Ri 2 (fc) both require 
k > 2. When fc = 2, Si (2) = (b 2 bi - bib 2 )Ds is in S?. 
When k = 3, Rn(3) = (bib 2 bi - b 2 Db?)Ds = R?[1,0], 
and R i2 (3) = bi[b 3 b 2 ]Ds = R?[0, 1], where b 3 >- b 2 . 



Consider Si (3). There are 3 elements led by variable b 2 : 
(b 2 bib 2 -biDb 2 )Ds, (b 2 bib 3 — bib 3 b 2 )Ds, and (b 2 b 3 bi — 
bib 2 b 3 )Ds. They all belong to S?. There are two other 



elements in Si (3): 
bib 3 b 2 )Ds. By 



(b 3 bib 2 — bib 2 b 3 )Ds and (b 3 b 2 bi 



b 3 b 2 bi 
b 3 bib 2 



bib 3 b 2 
bib 2 b 3 



G3 
G3 



bib 2 b 3 — b 2 b 3 bi; 

(b 2 b 3 bi — bib 2 b 3 ) 
+ (b 2 bib 3 — bib 3 b 2 ) 



both are in (S?) +1 D [M] 



So the statement in Step 1 holds for k < 3. Assume that 
it holds for all k < h. When fc = ft, we need to make I a - 
reduction to the leaders of the elements of any of the types 

S x (ft), Rii(ft), Ri 2 (ft), Ru [*](/»), Riat*](fc), (6.6) 

at the same time keep the reduction result to be within the 
K-linear space spanned by elements of the types listed in 
(|6.6p but where the left length h is replaced by all i < h. 
Then ultimately all the leaders of these elements become 
Z n -normal. 

Step 3. Consider types Rn(ft), Rn[*](/i), Ri 2 (ft), Ri2[*](ft). 
Let there be an Rn(ft)-typed element / = 2(bi[Abi]Ds), 
and an Ri 2 (/i)-typed element g = 2(bi [Bb 2 ]Ds), where B 
does not contain bi. In the following we omit the factor Ds. 

Do I D -reductions to A,B, and assume that the results are 

C' A b I); Da b E A b] A x . (( . ?) 



A 
B 



T- 



C B b 2 + b 2 D B + b 2 E B b 2 + B 



.v 



where (i) none of the terms in CajDajEa^" has bi at 
any end; 

(ii) none of the terms in Cb, Db, Eb, B^ has b 2 at any end; 

(iii) any of the four terms in each result may not occur; 

(iv) the component on the right side of the square symbol 
in each term, together with the symbol itself, are omitted, 
as they do not affect the analysis below; 

(v) in the extreme case, or B^ may be in K, if all 
vector variables in the term form squares and are moved to 
the right side of the square symbol. 

Substituting the reduction results into /, g, we get 
/ : 



(D A bi - (-l) fc D^bi)Db? 


G (Si(ft- 


-2), 




Ri(ft- 


-2)) 


+(biC A -(-i)' 1 biCj v )nb? 


G Ri(ft- 


2) 


+ (E A Db? - (-l) h biE A bi)Db? 


G Rii(ft- 


-2) 


+2(b 1 [A Jv bi]), 


G Rii(ft) 




2(bi[C B ])Db| 


G Ri(h- 


2) 


+2(bi[D B ]pb! 


G Ri(h- 


2) 


+2(b 1 [E B b 2 ])Db| 


G Ri 2 (ft- 


-2) 


+2(b 1 [B Jv b 2 ]). 


G Ri 2 (/i) 





9 = 



Notice that the left lengths indicated on the right column are 
the maximal possible ones for the corresponding types. So 



by induction hypothesis, we can assume that in /, g, mono- 
mials A = A^, B = B^ and both are I D -normal. 

The X n -reduction to the leaders of f,g axe much the same 
with the procedure in the proof of Theorem [7] starting from 
Step 3 there to Step 8, with negligible revisions. Formula 
(14.111) can also be used to split the leader of factor bi[Yibi 
■ ■ ■ Yfcbi] in a type-Rii[*] element, and the leader of factor 
bi[Yib 2 ■ ■ • Yfeb 2 ] in a type-Ri 2 [*] element. 

By induction on the order of the leader, we get that Rn(ft), 
Ri 2 (/i),Sq 1 (/i),Sq 1 [*](/i),Rn[*](/i), Ri 2 [*](ft) are all in 
where k — h. 

Step 4- Consider a general type-Sio(ft) element / = (Abi — 
biA)Ds. Let the I n -reduction result of A be as in (16. 7[) . 
Then if omitting "Lis", 



C A Db? - 
+biD A bi 
+ (biE A - 
+A Jv bi 



biC A bi 
- D A Db? 
-E A bi)Db? 



G Sqi(/i) 
G Sqi(ft) 
G Si(ft-2) 
G Sio(A) 



So we can assume that in / = (Abi — biA)Ds, monomial 
A = A^ and is I n -normal. 

The I D -reduction to the leading term of / is much the same 
with the procedure in the proof of Lemma [S] starting from 
Case (i) there to Case (ii). By induction on the order of the 
leading term, we get that Sio(fe) is in (|6.5p where k = h. 

Step 5. Consider a general type-Si(fe) element g = (AbiB — 
biBA)Ds, where A, B are both non-empty. Let the X n - 
reduction results of A,B be as in (|6.7p . where every b 2 is 
replaced by bi. Then if omitting "Ds", 

g = C A C B biDb?-biC B biC A bi 
+ (C A biD B - D B C A bi)Db? 
+ (C A biE B bi - E B biC A bi)Dbf 
+C A B JV Db? - biB^CA.bi 
+biD A biC B bi biC B D A Dbf 
+ (biD A D B D B biD A )Db? 
+biD A E B biDb? - E B D A Dbf 
+biD A biB iv biB Jv biD A 
+ (biE A C B bi - biC B E A bi)Db? 
+ (biE A biD B - D B biE A bi)Db? 
+biE A biE B biDbf - E B E A biDbf 
+biE A B Jv bf biB^biEA.bi 
-fA^biCBbi - biCebiA^ 
+ (A"D B - D B A w )Db? 
+A Jv biE B bi - biEebiA^ 
+A JV biB JV - bxB^A^. 

In the above result, the lines that do not belong to the ideal 
(Sq^ft), Si (ft - 2), S qi (ft - 2), Si (ft - 4)) are 



biD A biB JV biB^biDj^ 
+A w biC B bi - biCsbiA^ 
+A Jv biE B bi - biEsbiA^ 
+A Jv biB JV — biB^A^. 



G (Sii(ft),S qi (ft)> 
G (Sn(ft),Sio(ft)) 
G (Sn(ft),Sio(ft)) 
G Si(ft) 



(6.8) 



Some remarks on (|6.8|l are necessary. The first line of ( 



nonzero, is Sn(ft) when B^ ^ K, and Sq-^/i) otherwise. By 
A^biCebi-biCBbiA^ = (A Ar biC B bi-biA Jv b 1 C B ) + 
(bi A^biCe - biC B bi A^), the second line of (EH]) is a K- 
linear combination of an element of type Sio(/i) and another 
element of type Sn(/i). 

Consider a general type-Sn(ft) element p — (biAbiB — 
biBbiA)Ds, where the length of A is a > 0. When omitting 
□s, 

p T = biB(biA- (-l) a A t bi) 

+(-l) a A t BDb? - biBbjA 

= (-l) Q (A + BDbf - bjBA+bi) 

Sigh) ( _ 1)Q(A f B _ BA t) Db 2 

G (Si(ft-2)>. 
So Sn(ft) is in ()6.5p where k — h. 

By (|6.8|) . we can assume that in g = (AbiB — biBA)Ds, 
monomials A = A^, B = B^ and both are I D -normal. The 
I n -reduction to the leading term of g is much the same with 
that in the proof of Lemma [5] starting from Case (i) there 
to Case (iii). By induction on the order of the leading term, 

we get that Si(h) is in (|6.5[) where k — h. 

Consider a bracket polynomial / whose multiset of variables 
is A4. In any term of /, when all the brackets but one are 
expanded into two terms by definition, / is changed into a 
uni-bracket polynomial g. Using the Grobner base BG[M] 
to make reduction to g results in a uni-bracket polynomial h, 
where each term is 2/ -normal, h must have the lowest order 
lexicographically among all uni-bracket polynomials equal to 
/. It is called the lowest-representative normal form of g, or 
the uni-bracket normal form of /. 

Remark. In the above definition of normal forms, we only 
considered square-free ones. Of course any normal form can 
be converted to a canonical one, where all representatives 
are I-normal instead of Z D -normal. Later on, we consider 
only square-free ones. 

Let the size of M be m. Given any partition (ii, . . . , ijt) 
of integer m, where each ij > 1, there is a Gaianiello ex- 
pansion [5] of uni-bracket polynomials into bracket polyno- 
mials where each term is composed of k brackets of length 
ii, . . . ,ik respectively. Each expansion produces a normal 
form. It is not clear if such a normal form is of any value. 

7. NORMALIZATION OF BRACKET 
POLYNOMIALS 

In this section, we do NOT remove bracket symbols from 
uni-brackets in multiset of variables M. 

Consider attaching an additional vector variable vo ^ M to 
M to form a bigger multiset M. Let vo -< all variables in 
M . In the procedure of obtaining by the X n -reduction of 
bi[A] in the proof of Theorem 1 121 or in more details, in the 
proof of Theorem[?l if the input is vo [A] where A G Q m [M], 
then among the Grobner base of X n [M], only those elements 
in C/ D [7V(] are needed in I D -reduction. 



The reduction of vo[A] is a procedure of recursively do- 
ing I n -reductions to the leader of the bracket polynomial 
obtained from previous I n -reductions to [A]. At any in- 
stance, the representative of a bracket monomial in reduc- 
tion is the leader of the bracket monomial. The reduction 
results in a K-linear combination of monomials of the form 
vo[YiZi] • • ■ [YfcZfc], where z^Yi is ascending for every i. Af- 
ter removing vq from the result, we get another normal form 
of uni-bracket [A] in A4 . 

Theorem 13. Let f be a square-free bracket polynomial 
in multiset of variables M. Do the following to f: 

1. Always select the leader of a bracket as its representative. 

2. Rearrange the order of the brackets in the same term, so 
that the leading variables of the brackets are non- descending. 

3. Use \4-b]l to normalize the interior of a bracket; it also 
splits a bracket into two. 

4- Use ^4-9$ t o absorb a bracket into the one ahead of it. 

5. Use $4-10]} t° decrease the order of the product of two 
brackets by lifting a lower-order variable from the second 
bracket to the first. 

6. Use \4-H\l to segment a long bracket of type Y1Z1Y2Z2 ■ • ■ 
YfcZfc into short ones. 

7. Once the representative of the leading term of f is lP- 
normal, output the leading term, and continue the above T n - 
normalization to the remainder of f. 

The output, called the lowest-leader normal form, or leader- 
normal form is a bracket polynomial where the representative 
of each term is up to coefficient of the form [Y1Z1] [Y2Z2] • ■ ■ 
[YfcZfcJDs, where each z^Yi is ascending, Y1Y2 ■ • ■ Yfc and 
Z1Z2 • ■ • Zfc are both non-descending, and the leading variable 
U of each Y; satisfies L. >- Zj . Two bracket polynomials are 
equal if and only if their leader-normal forms are identical. 

Proof. After operations 1 and 2, for any bracket mono- 
mial in the reduction procedure, its representative is also 
its leader, so that the representative of the leading term of 
bracket polynomial / is the leader of /. Once the leader of 
/ is X D -normal, it is the leading term of the normal form 
of vector-variable polynomial / with respect to the Grob- 
ner base £/ n [A4]. By induction on the order of the output 
terms from the highest down, we get the uniqueness of the 
leader-normal form for /. 

For two equal bracket polynomials, they have identical uni- 
bracket normal forms, and so have identical leader-normal 

forms. 

From the above proof, we see that the leader-normal form 
of a bracket polynomial / has the following properties: (1) 
the representative of any term is the leader of the term; (2) 
the representative of the leading term is the leading term of 
the I D -normal form of vector-variable polynomial /. 

By Theorem[T2l the Grobner base BG[M] of J n [M] is com- 
posed of C/ n [A4] and the R?[*] : vagOs for all monomials 
g in the leader-normal forms of bracket polynomials in M, 
such that gs has length m. This phenomenon is easy to 
understand: the K-linear subspace of fT] D [.M] composed of 
polynomials whose terms are led by variable vo, is the space 



of degree-(m + 1) polynomials of the form Vo/Ds, for all 
bracket polynomials / in A4 such that fs has length m. 
The leader-normal forms are a basis of the K-linear space of 
length-m square- free bracket polynomials in M. 

In a leader- normal form, if we commute Y; and Zj in [YiZj], 
we get another normal form whose terms are up to coefficient 
of the form [ziYi] • • ■ [z^YfeJOs. If we write the left side 
of the square symbol in the following tableau form, where 
Yj = y l iyi 2 • ■ ■ y*t, , we get 



zi yn 

Z2 y2i 



yit! 



y2t 2 



ykt k 



where (1) each row does not need to have equal length, and 
it is not required that the length be non-increasing as in 
Young tableau; 

(2) each row is an ascending sequence of variables; 

(3) each column is a non-descending sequence of variables; 

(4) yu z < y(i+i)i for 1 < i < k. 

Such a normal form is called the straight form. Feature (4) 
above makes this definition stronger than the straight form 
(or standard form) of Young tableau. In comparison, in 
classical bracket algebra a bracket monomial is in straight 
form if and only if the entries are ascending along each row, 
and non-descending along each column. 

The procedure of deriving the straight form of a bracket 
polynomial is called straightening. Among the formulas used 
in Theorem [13] for straightening, (|4.10|l is highly nontrivial 
and requires further investigation. 



Set Bu, aD in (14.10P to be new A, B respectively, and let 
the lengths of B,C be b,c. Then (|4~K)]l can be written 
succinctly as follows: 

[Av][BwC] = [AwCv][B] - (-l) !l [AwB t v][C] 

-(-l)VCv][ABt] + (-l)VBV][AC] 
-[Aw][CvB]. 

(7.1) 

It is called the shuffle formula for bracket normalization. 

Proposition 14. For any two monomials A, B of length 
a, b respectively, 

AB + BA 



[AB] + (-l) 6 (A[Bt]-[A]Bt) 



AB — (— l) a+6 A t B t 



(-l)«([At]B-At[B]) 



(7.2) 



Proof. 



AB + BA 

2[AB] - (-l) a+6 B t A t + 2[B]A - (-l) 6 B f A 
2([AB] + [B]A-(-l) 6 B t [A]). 



□ 



The fundamental I-reduction formula (14.5[l is a direct con- 
sequence of the first identity in (|7.2|l for AB = uD. The 



shuffle formula (|7.1[) is a consequence of the following iden- 
tity by making left multiplication with A and then applying 
the bracket operator to both sides of the identity: 

v[BwC] + w[CvB] = (-^(wCvlB'] - [wCv]B' 

-wB f v[C] + [wB f v]C). 

(7.3) 

The identity can be obtained as follows: by the second iden- 
tity of (|7.2[) from right to left, the right side of (|7.3[) equals 



2' 1 ( BwC v + v BwC ) - (-l) 6+c 2~ 1 ( B t vC t w + w B'vC 1 ), 

which by the first equality of (|7.2I) . equals 

[vBwC] - (-l) i,+c v[C t wBt] 

- (-l) f,+c [wB t vC t ] +w[CvB] 
= v[BwC] + w[CvB]. 

The shuffle formula can be further generalized. In monomial 
[AvD][BwC], let the leading variable and trailing variable 
of any sequence F be If and tF respectively. Assume 1a ^ 
1b, and 1a >- tu, and 1b >- tc. Further assume 1a >- v >- 
w. Then X D -reductions can be made to [AvD][BwC] to 
decrease its leader, leading to the following result: 



Proposition 15. Let the lengths of monomials A, B, C, D 
be a, b, c, d respectively. Then 

[AvD][BwC] = [vDBw][AC] - (-l) Hc [vDC t w][AB t ] 
-(-l) d [Aw][D t vBC] - [vD][AwBC] 
-(-l) 6 [AwB f vD][C] + [AwCvD][B]. 

(7.4) 



Proof. 

4[AvD][BwC] 

= A vDBw C (-!) ! ' +c A vDClw B t 

-(-l) a+d D t vA t BwC + (-l) a+6+c+d D t vA t C t wB t 

X = (vDBw+ (-l) 6+d wB t D t v)AC 

-(-l) b ((-l) c vDC + w + (-l) d wCD t v)AB t 
_(_ 1 )6+d AwB t D t vC + (_i)Hi AwCD t vB t 
_ ( _ 1) a+d D t v(BwC _ (_!)Hc c t wB t) A t 

Z = 4([vDBw][AC] - (-l) f,+ °[vDC t w][AB t ]) 
+(-l) a+d (BCD t v - D t vCB)wA t 
+(-l) HJ Aw(CD t vB t - BtD+vC). 

(7.5) 

By (EH), 

(-l) d (BCD t v - D+vCB) 

= (— l) d BCD t v - (-l^vDCtBt 
+(-l) b+c (vD - (-l) d D t v)C t B t 
+(-l) c+d 2 {D t v(-[C t ]B + C+fB])}, 

we get from (|7.5p the following: 

4[AvD][BwC] 

z = 4([vDBw][AC] - (-^^[vDC+wlIABt]) 
+4 [{(-l) d [BCD t v] + (-l) i,+c [vD]C t B t 
-(-l) c+d D t vB[C t ] + (-l) c+d D t vC t [B]}(-l) a wA t ] 
+(-l) b Aw{-(-l) c vDC t B t + (-^"BtC+vD 
+(-l) d CD t v B t - (-l) d B t D t vC} 



4([vDBw][AC] - (-l^+^vDC+wHABt]) 
-4{(-l) d [Aw][BCD t v] + [AwBC][vD] 

+(-l)'[AwB t vD][C] - [AwCvD][B]}. 



+ 



□ 



Clearly each term in the result of (|7.4[) has lower leader than 
the input. To better understand this reduction formula, we 
write it in tableau form: 



A((-l) 6+1 Bw) t vD 
C 

vD 

((-i) 6 Bt)w((-irA)t ( (-irct) . 

wC((-l) d+1 vD)t 
B((-l) a A)t 



AvD " 




" AwCvD 


BwC 




B 



+ 



+ 



+ 



((-l) i,+1 Bw) t ((-l) d+1 vD) t 
AC 



+ 



Aw 

((-l) b Bt)vD((-l) c Ct) 



(7.6) 



As 1a >~ v y w and 1a >- tu, to decrease the order of the 
leader, in the first line of (|7.6[1 . a subsequence or reversed 
subsequence of the second row of the input tableau is moved 
up between A, v of the first row, with the requirement that 
the subsequence be led by w. In the second line of (1T.6|) . A 
is moved down behind w of the second row. 



In the third and fourth lines of (|7.6|) . the leading subsequence 
A of the first row is commuted with a subsequence or re- 
versed subsequence of the second row led by w. In the last 
line of (|7.6[1 . the trailing subsequence vD of the first row is 
commuted with variable w of the second row. 

We make comparison with the shuffle formula for straighten- 
ing in classical bracket algebra [13] . In the classical bracket 
algebra over 3D vector space, where the exterior product is 
also denoted by juxtaposition of elements, suppose that 



a v d 
b w c 



is not straight: avd is ascending, so is bwc; ab is non- 
descending, but v >- w. The shuffle formula is obtained 
as follows. For any four vectors v, d, b, w of the 3D vector 
space, their exterior product equals zero. By 



= a V (vdbw) V c 

= [avd] [bwc] — [avb][dwc] + [avw][dbc] 
+ [adb][vwc] — [adw][vbc] + [abw][vdc] 



(7.7) 



where "V" is the dual of the exterior product called the meet 
product |14] . we get the following shuffle formula, also called 
van der Waerden relation: 



a 


V 


d 




a 


V 


b " 




a 


v 


w 


b 


w 


c 




d 


w 


c 


+ 


b 


d 


c 



a 


b 


d 


+ 


a 


w 


d " 


V 


w 


c 


b 


V 


c 



a b w 



(7.8) 



(1) The first line commutes d of the first row and one of the 
first two vectors of the second row; 

(2) the second line commutes v of the first row and one of 
the first two vectors of the second row; 

(3) the last line commutes vd of the first row and the first 
two vectors of the second row. 

For straightening in classical bracket algebra, there are other 
operations besides (I7.8[) . When a y b, we only need to 
commute the two brackets. When d y c while a ^< b and 
v ■< w, we have another formula for straightening. By 

= (av) V (dbwc) 

= [avd] [bwc] — [avb] [dwc] 
+ [avw][dbc] — [avc][dbw], 

we get the following Grassmann-Pliicker relation: 





a 


V 


d 




a 


V 


b " 




a 


V 


w 




a 


V 


c 




b 


w 


c 




d 


w 


c 


+ 


b 


d 


c 


+ 


b 


w 


d 



The last vector d of the first row commutes in turn with 
every vector of the second row. 

8. CONCLUSION 

In the bottom-up approach to manipulating brackets, long 
brackets are expanded into basic ones by Caianiello expan- 
sion, in the end only brackets of length 2 and 3 are left 
for further algebraic manipulations. This approach proves 
to be inefficient in practice, despite the fact that there are 
straightening algorithms for polynomials of basic invariants 

Uni-bracket polynomials provide a top-down approach to 
manipulating brackets. Given a bracket polynomial, by "un- 
grading", each bracket but one in every term is expanded 
into a vector-variable binomial, and the bracket polynomial 
is changed into a uni-bracket one. Algebraic manipulations 
of uni-bracket polynomials can take full advantage of the 
associativity of the vector- variable product and the symme- 
tries within a uni-bracket. The Grobner base C/ D [A1] pro- 
vided by this paper further fulfills the arsenal of symbolic 
manipulations on uni-bracket polynomials. 

The last section of this paper suggests a third approach to 
manipulating brackets by algebraic manipulations directly 
upon the input brackets. To establish this approach there 
are many research topics ahead: division among bracket 
polynomials, properties of principal ideals, bracket polyno- 
mial factorization, and simplification by reducing the num- 
ber of terms, etc. This seems to be a promising approach. 

This paper is supported partially by NSFC 10871195, 
60821002/F02, and NCMIS of CAS. The Grobner bases for 
m = 5, 6 were first computed by Dr. L. Huang on his imple- 
mentation of the classical non-commutative Grobner base 
algorithm [9]. 
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